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Abstract 

Absorption of radiation from the sodium dimer molecular states correlating 
to Na(3s)-Na(3s) is investigated theoretically. Vibrational bound and contin- 
uum transitions from the singlet X 1 state to the first excited A 1 and 
B 1 !^ states and from the triplet a 3 £+ state to the first excited b 3 X+ and 
c 3 IT 9 states are studied quantum-mechanically. Theoretical and experimen- 
tal data are used to characterize the molecular properties taking advantage 
of knowledge recently obtained from ab initio calculations, spectroscopy, and 
ultra-cold atom collision studies. The quantum-mechanical calculations are 
carried out for temperatures in the range from 500 to 3000 K and are com- 
pared with previous calculations and measurements where available. 
PACS numbers: 33.20.-t, 34.20.Mq, 52.25.Qt 
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I. INTRODUCTION 



Vast amounts of experimental spectroscopic data on the electronic states and ro- 
vibrational levels of the sodium dimer are available and many theoretical studies have been 
performed. For example, Ref. |1J presents an extensive bibliography summarizing a variety 
of work dating from 1874 to 1983. Nevertheless, recent developments in atom trapping and 
cold atom spectroscopy have led to improved atomic and molecular data through combi- 
nations of cold collision data, photoassociation spectroscopy, and magnetic field induced 
Feshbach resonance data [0-0. 

Now that very reliable information is available, calculations of absorption spectra at high 
temperatures become feasible. Absorption coefficients in absolute units for a gas of sodium 
atoms and molecules at temperatures from 1070 to 1470 K were measured over the range 
of wavelengths from 350 to 1075 nm by Schlejen et al. ||. They performed semiclassical 
calculations involving the relevant molecular singlet and triplet transitions, however, those 
previous calculations do not fully reproduce their experimental spectra 0. The present 
work is concerned with the absorption involving two ground Na (3s) atoms and a ground 
Na (3s) atom and an excited Na (3p) atom, corresponding to transitions between the singlet 
transitions from the X 1 S+ state to the A X S+ and B 1 U U states and the triplet transitions 
from the a 3 S+ to the b 3 S+ and c 3 IT g states. We assembled and evaluated the available 
data for the molecular system and calculated quant um-mechanically the absorption spectra 
at temperatures between 500 and 3000 K. 

II. ABSORPTION COEFFICIENTS 

The thermally averaged absorption coefficients k v for molecular spectra at wavelength v 
are obtained from the product of the thermally averaged cross sections and the molecular 
density [T(J. In turn, the molecular density can be expressed in terms of the atomic density 
squared and the chemical equilibrium constant JTTJ . In the present study, we use the atomic 
density- independent reduced absorption coefficient, ||, where n a is atomic density. 

Four possible types of vibrational transitions between two electronic states can be iden- 
tified: bound-bound (bb), bound-free (bf), free-bound (fb) and free-free (ff) and quantum- 
mechanical expressions for the reduced absorption coefficient can be derived. The radial 
wave function for a bound level, with vibrational and rotational quantum numbers, v and 
J, is obtained from the Schrodinger equation for the relative motion of the nuclei, 
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d 2 <P v jA(R) , ( n jp J(J+1)-A 2 \ 

— — + I 2^E vJA - 2fiV{R) — 2 I (pvJA{R) = 0, (1) 

where V(R) is the potential for the relevant electronic state labeled by the projection A of 
the electronic orbital angular momentum on the internuclear axis, /x is the reduced mass of 
the nuclei, and E is the eigenvalue of the bound level or the continuum energy. 

For the temperatures of interest here, T < 3000 K, the bound-bound reduced absorption 
coefficient from a vibration-rotation state of the lower electronic state (v", J", A") to the 



vibration-rotation state of the upper electronic state (v', J', A') is p!0| . |T2| 

% = ^T-I^bT) exp(D e /k B T) 
n 2 a h 

x E E^'^' + ^exPt-^W^)!^ (2) 

v" J" v' 

where v is the frequency, v is the transition energy of the bound-bound transition, 

n ,s. _ (2 - a , A '+A") 87T 3 ?y 

C[V) ~ 2 - 5 0A „ 3^ (3) 
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S m and S a are spin multiplicities for, respectively, the Na molecule and the Na atom, and ks 
is Boltzmann constant. The Q-branch approximation (J' = J") is used and the line-shape 
function g[y — v) is replaced by 1/Au. In evaluating Eq. (Q) at some Ui on the discretized 
frequency interval, all transitions within the frequency range Ui — \Av to Vi + |Az/ are 
summed to give a value k ^ . The nuclear spin statistical factor ojj for 7 Na 2 with I = | is 
[l/{21 + 1)] = | for even J and [(J + 1)/(2J + 1)] = § for odd J. 

The bound-free absorption coefficient from a bound level of the lower electronic state 
(v"J"A") to a continuum level of the upper electronic state (e'J'A') is 

k h j 



^ = C(v)f(k B T)exp(D e /k B T) 



x J2 UJ"(2J" + 1) exp(-Ev»j»/k B T) \{4>v"j"K"\ D(R) |0 e vw)| • (5) 

v" J" 

The continuum wave function is energy normalized. For a free-bound transition and free-free 
transition, respectively, 

-\ = C{v)f{k B T) 



n 2 



x Uj»(2J" + l)exp(-e"/k B T) \(4> e "j"A"\ D(R) |^j'A')| 2 (6) 

v' J 1 
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III. MOLECULAR DATA 



The adopted singlet X 1 S+, A 1 S+, and B 1 TL U potentials and the differences of the upper 
potential and the lower X 1 E+ potential (difference potentials or transition energies) are 
plotted in Fig. [I]. The adopted triplet a 3 S+, b 3 E+, and c 3 IT 9 potentials and the difference 
potentials are plotted in Fig. ^|. In the remainder of the section details on the construction 
of the potentials are given. We use atomic units throughout. 



A. The X : S+ potential 

For R < 4 ao, we adopted a short range form aexp(— bR), with a = 2 702 514.0 cm -1 and 
b = 2.797 131 A -1 as given by Zemke and Stwalley [|16j. Over the range of 4 < R < 30 ao we 
used the Inverse Perturbation Analysis (IPA) potential given by van Abeelen and Verhaar J7| 
which is consistent with data from photoassociation spectroscopy, molecular spectroscopy, 
and magnetic-field induced Feshbach resonances in ultra-cold atom collisions. For the long 
range form, we used 

- C 6 /R 6 - C 8 /R 8 - C 10 /R 10 - AR^- 1 exp(-2a J R), (8) 

where C 6 = 1561 @, C 8 = 111877, C 10 = 11065 000 |7|, A = and a = 0.626 @JH. To 
fit the very accurate dissociation energy, 6 022.0286(53) cm -1 , recently measured by Jones et 



al. ||19|| , a point at the potential minimum 5.819 460 ao was added. The short and long-range 
data were smoothly connected to the IPA values. Vibrational eigenvalues calculated with 
our adopted potential agree for v < 44 to within 0.1 cm -1 with published Rydberg-Klein- 
Rees (RKR) values ||16|| . Our final potential yielded an s-wave scattering length of 15 ao in 



satisfactory agreement with the accepted value of 19.1 ± 2.1 ao 



B. The A X E+ potential 

We used ab initio calculations given by Konowalow et al. [OT for values of R over the 
range 3.8 a < R < 4.75 a . We combined the RKR potential values over the range 
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2.522 19 A< R < 7.204 14 A given by Gerber and Moller pTJ with the RKR potential 
values over the ranee of 7.260 536 A< R < 261.327403 A given by Tiemann, Knockel, and 



Richling 22,23 . The data was connected to the long range form, 



— C 3 / R 3 — CqJ R 6 — C 8 / R 8 



(9) 



with the values of C 3 = 12.26, C 6 = 4 094 and C 8 = 702 500 ||]. For R < 3.8 a , the form 
aexp(— bR) + c was used with the parameters a = 0.9532, b = 0.5061 and c = 0.104696 
computed to smoothly connect to the RKR points. The adopted potential yields a value of 
D e = 8 297.5 cm -1 using T e = 14 680.682 cm -1 and the atomic asymptotic energy of 
16 956.172 cm -1 [25]. The calculated eigenvalues reproduce the input RKR values to within 
0.4 cm -1 . For the transition frequencies measured by Verma, Vu, and Stwalley [|26| and by 
Verma et al. over a range of vibrational bands we find typical agreement to about 0.4 cm" 1 
for J' values up to 50 increasing to 1 cm -1 for J' = 87. We also have good agreement with 
less accurate measurements by Itoh et al. [|27]]. One precise transition energy measurement 
is available: In a determination of the dissociation energy of the sodium molecule Jones et 
al. 0] measured the value 18762.3902(30) cm" 1 for the v' = 165, J' = 1 to v" = 31, J" = 
transition energy. Our value of 18762.372 cm -1 is in excellent agreement. 



C. The B l U u potential 



The RKR potential of Kusch and Hessel [28] was used[] over the range of 2.655 5671 A < 
R < 5.173 513 4 A. For the values of R in the ranges 2.581 A < R < 2.646 0268 A and 



5.251 918 4 A < R < 11.0 A, we took the potential values from Tiemann |30| . We also took 
his long-range form, 



C 3 /R 3 - C 6 /R 6 + C 8 /R 8 - aexp(-bR), 



(10) 



with C 3 = 6.1486, C 6 = 6490.5, C 8 = 868135.2, a = 23.7011, and b = 0.7885. For 
R < 2.581 A, the form aexp(-bR) + c was used with the values a = 14.97332, b = 1.42983 
and c = 0.0121935 chosen to give a smooth connection with the data from Tiemann. 

The B l U u potential exhibits a barrier that has been studied extensively PT| , P9"| -|3"T[] and 
the maximum value occurs around R = 13 ao (6.9 A) as shown in Fig. [l|. We took D e = 



For this reference, we correct an apparent typographical error of 4.309 78 A with 4.339 78 A ob- 
tained by comparison with RKR potential of Demtroder and Stock |E9]. 
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2 676.16 cm 1 using T e = 20 319.19 cm 1 from Kusch and Hessel [p8 and the barrier energy 



371.93 cm 1 measured from dissociation given by Tiemann pD |. The calculated energy 



23 393.524 cm -1 of the v' = 31, J' = 42 state with respect to the X 1 S+ state potential 
minimum compares well to the measured value, 23 393.650 cm' 1 . Quasibound levels from 
v' = 24 to v' = 33 for the several J' values observed by Vedder et al. are reproduced 
to within 0.1 cm -1 and calculated transition frequencies compare well, to within 0.5 cm -1 , 
with those measured by Camacho et al. 



D. The a 3 £+ potential 
RKR potentials are available from Li, Rice, and Field MM and Friedman-Hill and 



Field [35] and a hybrid potential was constructed by Zemke and Stwalley [16] using var- 
ious available data. An accurate ab initio study was carried out by Gutowski [[J for R 
values between 2 and 12.1 A and the resulting potential has well depth 176.173 cm -1 and 
equilibrium distance 5.204 A. 

Our adopted potential consists of Gutowski's potential connected to the long-range form 
given in Eq. (|8[) with the values for Cq, Cg, C\o and a the same as for the X X E+ state, but 
with A = — g^j. For R < 2 A the short range form aexp(— bR) was used with a = 1.4956 
and b = 0.79438 chosen to smoothly connect to the ab initio data. Our adopted potential 
yields an s-wave scattering length of 65 a$ in agreement with the value 65.3 ± 0.9 of van 
Abeleen and Verhaar 0. Recently, a potential alternative to Gutowski's was presented by 



Ho et al. |3"6|. For the present study the two potentials are comparable — we will explore 



their differences in a subsequent publication. 

E. The b 3 S+ and c 3 n 9 potentials 

We are unaware of empirical excited state triplet potentials but ab initio calculations 
are available from Magnier et al. f3"7[| , Jeung and Konowalow et al. ||2U|| . Comparing the 
available potentials, we found for the b 3 £+ state that the experimentally determined 
T e of 18 240.5 cm -1 and D e of 4 755 cm -1 are closest to Magnier's calculated values (T e of 
18 117 cm -1 and D e of 4 740.7 cm -1 ) compared with Jeung's (T e of 18 400 cm -1 and D e of 
4 702.4 cm -1 ) and Konowalow et a/.'s (D e of 4 599 cm -1 ). Also, we found Magnier's potential 



gave the best agreement with experimental measurements |4(| of the term differences of the 
a 3 E+(V') — > b 3 £+(V) vibrational transitions. For the b 3 S+ state and, in the absence of 
experimental data for the c 3 n„ potential, we adopted Magnier's calculated potentials over 
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the range of R values 5 < R < 52 a . Over the range of R values 4.25 < R < 5 ao, we used 
potentials by Konowalow et al. |[20|| . For the b 3 E+ and c 3 Il 9 adopted potentials, the long- 
range form was taken from Marinescu and Dalgarno |24j for R > 52 a and for R < 4.25 ao, 
we used the form aexp(— bR) where the values are a = 55.7864 and b = 1.75934 for the 
b 3 S+ potential and a = 2.67691 and b = 0.91547 for the c 3 IT 9 potential. 



F. Transition dipole moment functions 



We used for the singlet transitions the ab initio calculations of Stevens et al. [[|l[ over the 
range 2 < R < 12 a . For R > 12 the transition dipole moment functions were approximated 
by a + b/r 3 , where a = 3.586 4 and b = 284.26 for X X E+ — > A 1 ^^ transitions and a = 3.5017 
and b = —142.13 for X x £+ — > B 1 U U transitions. The parameter values for a were selected 
to match the short-range parts and those for b were from Marinescu and Dalgarno []24|. The 



X — > A dipole moment function was scaled with a factor of 1.008, as discussed 



in Sec. [IV A| below. For the triplet transitions the ab initio calculations of Konowalow et 



al. [42] were used over the range 4 < R < 100 ao- 



IV. RESULTS 
A. Lifetimes 



In order to evaluate our assembled potential energy and transition dipole moment data 
we calculated lifetimes of ro-vibrational levels of the A X S+ and B 1 Il tJ states and compare 
with prior studies. 

Lifetimes for levels of the A X S+ state have been measured p6| , ^3| -|45| and calcu- 
lated ^6j,^6j. In Fig. H] we present a comparison of rotationally resolved lifetimes for levels 



of the A X S^ state measured by Baumgartner et al. [45] with the present calculations. In 
evaluating the lifetimes, we used the procedures described in Ref . [|H]] . When the transition 
dipole moment function of Ref. |41[] is multiplied by a factor of 1.008, agreement is generally 
very good over the range to 3500 cm -1 of available term energies. Our calculations are 



also in good agreement with the rotationally unresolved measurement of Ducas et al. |43| 



and the calculations using different molecular data by Pardo [46 



Rotationally resolved lifetimes for the B 1 I1 M state have been measured by Demtroder 
et al. P7| . Demtroder et al. found that the lifetimes for the B 1 I1 U state are sensitive to 
the slope of the transition dipole moment function in the range of internuclear separation 
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from, roughly, 4 < R < 10 ao and they obtained an empirical value for the function that 
we found to be in good agreement with the transition dipole moment function of Stevens et 
al. ||41[|. Using the transition dipole moment function of Stevens et al, in turn, we find good 



agreement between our calculated lifetimes and experimental lifetime measurements g7| , as 



shown in Fig. £|. The ab initio dipole moment of Konowalow et al. |42| was found not to 
reproduce the experimental lifetimes. 



B. Absorption Coefficients 

Absorption spectra in the far-line wings of the Na(3s)-Na(3p) resonance lines are investi- 
gated in terms of singlet and triplet molecular transitions. The blue wing consists of radiation 
from X X E+ — > B 1 U U and a 3 E+ — > c 3 II 9 transitions and the red wing from X — > A 1 E„ 
and a 3 £+ — > b 3 X+ transitions. There are few experimental studies P, |T4| , pE8"| | and that of 



Schlejen et al. is most relevant to our work. In this section, we compare our calculated 
absorption coefficients with the measurements of Schlejen et al. ||. 

The theoretical spectra are assembled from four molecular band spectra over the wave- 
length range 450-1000 nm excluding the region 589±2 nm around the atomic resonance 
lines. The far line wings are calculated using Eqs. (f2|)-(|7]), with the data from Sec. III. In 
the calculations, all the vibrational levels including quasi-bound levels with rotational quan- 
tum numbers up to 250 are included. The maximum internuclear distance that is used for 
integration of the transition dipole matrix element is approximately 100 ao and the Numerov 
integration used to obtain the energy-normalized continuum wave function is carried out to 
100 ao at which the wave function is matched to its asymptotic form. The bin size Av used 
for Eq. (|2|) was 10 cm -1 simulating the experimental resolution. Results for absolute absorp- 
tion coefficients computed with the quoted atomic densities and temperatures of Schlejen et 
al and shown in Fig. [5] compare very well with the four experimental spectra given by 
Schlejen et al, given in Figure. 5 of Ref. ||. The spectra show clearly that as temperature 
increases, certain satellite features grow more apparent at 551.5 nm and 804 nm. These 
satellites will be discussed in greater detail later in this section. The present calculations 
reproduce fine-scale ro-vibrational features present but unresolved in the measurements of 
Ref. §. 

We also have calculated reduced absorption coefficients at temperatures up to 3000 K 
using the bin size Au of 1 cm -1 for Eq. The contributions of the four molecular bands 
to the reduced absorption coefficients are shown in Fig. |6| for three temperatures 1000 K, 
2000 K and 3000 K. As can be seen by comparing columns (a) and (b) in Fig. || the singlet 



S 



transitions contribute more to the reduced absorption coefficients in the far line wings and 
the triplet transitions contribute more near the atomic resonance lines. We found that for 
singlet transitions bound-bound and bound-free transitions are dominant over free-bound 
and free-free transitions for the temperature range T < 3000 K, thus accounting for the 
"grassy" structure in Fig. 0(a). However, the free-bound and free-free contributions increase 
rapidly with temperature. In contrast to the singlet transitions, the triplet transitions arise 
mainly from free-bound and free-free transitions due to the shallow well of the initial a 3 S+ 
state. Hence, the reduced absorption coefficients in Fig. ||(b) do not exhibit much structure. 
Because the density of bound molecules decreases rapidly with increasing temperature, the 
reduced absorption coefficient in the line wings due to the singlet transitions also decreases 
rapidly with increasing temperature. It should be noted that the scale of the reduced 
absorption coefficient at 1000 K is two orders of magnitude larger than the scale shown for 
T = 2000 K and 3000 K. 



Woerdman and De Groot [48| derived the reduced absorption coefficient at 2000 K from 
a discharge spectra. The measured values of 5 ± 1 x 10 -37 cm" 1 at 500 nm and 10 ± 1 x 
10 -37 cm -1 at 551.5 nm are well-reproduced by our values of, respectively, 5 x 10~ 37 cm -1 
and 11 x 10~ 37 cm -1 calculated with the bin size Az/ of 5 cm -1 simulating the experimental 



resolution obtained by Woerdman and De Groot [48 



The molecular absorption spectra contain "satellite" features around the energies where 
the difference potentials possess local extrema [fl"5|J)C| ]. For Na 2 the satellite frequencies have 
been studied [^,|T^,^,^1|] and the energies have been calculated using ab initio methods E2 . 
In the present work, we investigate the satellites arising from a 3 S+ — > c 3 n g , X X E+ — > A 
and a 3 £+ — > b 3 £+ transitions with measured maximum intensities at, respectively, the 
wavelengths 551.5 nm, 804 nm, and 880 nm. The calculated extrema of the difference 
potentials adopted in the present study occur at wavelengths at 548 nm, 809 nm and 913 nm, 
however in the quantum-mechanical approach there is no well-defined singularity. 

We can use the quantum-mechanical theory to study satellite features in more detail and 
as a function of temperature. In Fig. |7| we show calculated reduced absorption coefficients 
at three temperatures for the a 3 £+ — > c 3 n 9 and X : S+ — > A X S+ transitions. The rich 
ro- vibrational structure in the X X S+ — > A 1 Ej" satellite feature arises because the dominant 

y a 

contributions are from bound-bound transitions; the structure is not reproduced by semi- 
classical theories 0. In contrast, the smooth, structureless a 3 S+ — * c 3 n 9 satellite feature is 
due mainly to free-free transitions, and consequently, the decrease of the satellite intensity 
with temperature is less severe. The slight discrepancy between the calculated wavelength 
of 550 nm and the measured wavelength of 551.5 nm [j48,52,53|l for the peak intensity is 
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probably due to remaining uncertainties in the triplet potentials |36|j54 . 

We also investigated the a 3 E+ — > b 3 £+ satellite which is far weaker in intensity at 
T < 3000 K than the X : S+ A*£+ and a 3 S+ -> c 3 n s satellites. The a 3 S+ -> b 3 E+ 
satellite arises primarily from free-bound transitions. The population density of atom pairs 
with high continuum energies in the initial a 3 S+ state increases with temperature, see 
Eq. (^|), and more ro- vibrational levels in the b 3 £+ state are accessible through absorption 
of radiation, as can be seen from the potential curves shown in Fig. 0(a). As a result, this 
satellite feature exhibits an increase in intensity with temperature. In Fig. [8] calculated 
absorption coefficients for temperatures 1000, 1500, 2000 and 3000 K are plotted. The 
satellite feature intensity was measured at 1470 K by Schlejen et al. [|]. They observed a 
primary peak at 880 nm and a secondary peak at 850 nm, compared to our calculated values 
at 1500 K of 890 nm and 860 nm, respectively. The 10 nm discrepancy in both peaks is 
probably a result of uncertainties in the short range parts of our adopted a 3 £+ and b 3 S+ 
potentials. Our calculations also demonstrate that the wavelengths of the peaks change with 
temperature, see Fig. and that the primary peak from quantum-mechanical calculations is 
less prominent than that obtained from semiclassical calculations exhibited in Figures 6(c) 
and 6(d) of Schlejen et al. ||. Our calculated reduced absorption coefficients appear to be 
in excellent agreement with the reduced absorption coefficients interpolated from Figures 
6(a) and 6(b) of Schlejen et al. using their quoted Na densities. 



V. CONCLUSIONS 



We have carried out quantum-mechanical calculations of the reduced absorption coeffi- 
cients in sodium vapor at high temperatures. Accurate molecular data are an important 
ingredient. Comparisons with experiments 0,|48| are good, but the theory is not limited 



by the previous experimental resolution. Future work |55| will focus on comparisons of the 
present theory and experiments currently on-going in our group 
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FIG. 1. (a) Adopted potentials V(R) for the X X S+, A 1 !!*, and B 1 !^ electronic states, (b) 
Difference potentials Vg ln (R) — V^_ 1 ^+(R) and Vj^ 1 - E +(R) — V^ 1 - E +(R). 
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FIG. 2. (a) Adopted potentials V(R) for the a 3 S+, b 3 E+ and c 3 II 9 states, (b) Difference 
potentials V h3j: +{R) - V a3j: +(R) and V C 3 Ug (R) - V & ^+{R). 
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(a) singlet bands 



(b) triplet bands 



(c) all bands 
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at 2000 K and 3000 K. The calculations were performed with bin size Au = 1 cm -1 . 
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